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Summary. — Some models of modified gravity and their observational manifesta- 
tions are considered. It is shown, that gravitating systems with mass density rising 
with time evolve to a singular state with infinite curvature scalar. The universe 
evolution durung the radiation dominated epoch is studied in _R^-extended gravity 
theory. Particle production rate by the oscillating curvature is calculated. Possible 
implications of the model for cosmological creation of non-thermal dark matter are 
discussed. 
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1. Introduction 

Discovering of the cosmic antigravity based on the accumulated astronomical data, 
such as observation of the large scale structure of the universe, measurements of the 
angular fluctuations of the cosmic microwave background radiation, determination of the 
universe age (for a review see [1]), and especially discovery of the dimming of distant 
Supernovae [2], is the most attractive event in cosmology of the last quarter of century. 
It was established and unambiguously proved that the universe expansion is accelerated, 
but the driving force behind this accelerated expansion in still unknown. 

Among possible explanations, the most popular is probably the assumption of a new 
(unknown) form of cosmological energy density with large negative pressure, P < — p/3, 
the so-called dark energy, for a review see e.g. [3]. 

A competing mechanism to describe the accelerated expansion is represented by grav- 
ity modifications at small curvature, the so-called /(i?)-gravity theories, as suggested in 
ref. [4] . In these theories the standard Einstein-Hilbert Lagrangian density, proportional 
to the scalar curvature R, is replaced by a function f{R), so the usual action of General 
Relativity acquires an additional term: 
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d^Xy^ [R + F{R)] + Sr, 



where nipi = 1.22 • 10 GeV is the Planck mass and Sm is the matter action. 
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The original version of such models [4] suffers from a strong instability in presence of 
gravitating bodies [5] and because of that more complicated functions F(R) have been 
proposed [6, 7, 8, 9], which are free from the mentioned exponential instability. 

Though free of instability [5], the models proposed in [6, 7, 8] possess another trou- 
blesome feature, namely in a cosmological situation they should evolve from a singular 
state in the past [10]. Moreover, it was found in refs. [11, 12] that in presence of matter, 
a singularity may arise in the future if the matter density rises with time; such future 
singularity is unavoidable, regardless of the initial conditions, and is reached in a time 
which is much shorter than the cosmological one. 



2. Explosive phenomena in modified gravity 

In paper [12] the model of modified gravity with F{R) function suggested in ref. [6] 
was considered: 



(2) 



F{R) = XRo 



R^ 



Here constant A is chosen to be positive to produce an accelerated cosmological expansion, 
n is a positive integer, and Rq is a constant with dimension of the curvature scalar. The 
latter is assumed to be of the order of the present day average curvature of the universe, 
i.e. Rq ^ ^/tfj, where tjy sa 4 • 10^^ sec is the universe age. 
The corresponding equations of motion have the form 



(m) 



(3) (1 + F') R^, -^(R + F) g^, + [g^^D^D^ - D^D,) F' = 



'Pi 



where F' = dF/dR, is the covariant derivative, and T^^f^ is the energy-momentum 
tensor of matter. 

By taking trace over ^ and v in eq. (3) we obtain the equation of motion which 
contains only the curvature scalar R and the trace of the energy-momentum tensor of 
matter: 



(4) 



iD'^F' - R + RF' - 2F = T . 



where T ^ 8nTi^/m%i. 

We analize evolution of R in massive objects with time varying mass density , pm 3> 
Pc- The cosmological energy density at the present time is pc ~ 10~^^g/cm^, while 
matter density of, say, a dust cloud in a galaxy could be about p„i ^ lO^^^g/cm'^. 
Since the magnitude of the curvature scalar is proportional to the mass density of a 
nonrelativistic system, we find R^ Rq. In this limit: 



(5) 



F{R) « -XRo 



1 - 



Ro 



2n 



(6) 



The equation of motion is very much simplified if we introduce the new notation 

n \ 2n+l 



w = -F' = 2nX ^ 
R 




Fig. 1. - Potential U{z) = z(l + kt) - ~ u), = i, r = 0. 



Evolution of w is governed by a simple equation of unharmonic oscillator: 

(7) (a^- A)w + c/'H = 0. 

Potential U (w) is equal to: 



(8) U{w)^^iT^2XRo)w + ^ 



where = l/(2n+ 1), q = 2nA, and in eq. (7) U'{w) = dU/dw. 

Notice that infinite R corresponds to w = —F' = 0, so if F' reaches zero, it would 
mean that R becomes infinitly large. 

Potential U would depend upon time, if the mass density of the object under scrutiny 
changes with time, T = T{t). If only the dominant terms are retained and if the space 
derivatives are neglected, equations (7), (8) simplifie to: 

(9) z" - + (1 + kt) = . 

Here we introduced dimensionless quantities 



Minimum of the potential U{z) (Fig. 1) sits at Zmin = (1 + kt)^^/"^. When the mass 
density rises, the minimum moves towards zero and becomes less deep. If at the process 
of "lifting" of the potential z(t) happens to be at ?7 > it would overjump potential 
which is equal to zero at z = 0. In other words, z(t) would reach zero, which corresponds 
to infinite R, and so the singularity can be reached in finite time (see Fig. 2). 

The simplest way to avoid singularity is to introduce i?^-tcrm into the gravitational 
action: 

(11) SF{R) = -R^/6m^, 



where to is a constant parameter with dimension of mass. 
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Fig. 2. - Ratio z{t) / z,nin{T) (left) and functions z(r) and Zmin{T) (right) for n — 2, n 
pm/pc = 10*. The initial conditions: z;(0) = 1 and z'{0) = 0. 
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In the homogeneous case and in the limit of large ratio R/Rq equation of motion for 
R is modified as: 
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R 



2n+2 



6\n{2n+l)R^ 



2n^ 



R - {2n + 2) 



fl^"+^(i? + r) 
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R 6Xn{2n + l)i?f 



With dimcnsionless curvature and time 



(13) 



R 



T, 
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6An(2n+ l)i?;^ 



1/2 



equation for R is transformed into: 

(14) (1 + .gz/"+2) y" _ 2(n + 1) + [y-{l + k.t,)] = , 

y 

where prime means derivative with respect to ti. 

Wc introduced here the new parameter, g, which can prevent from the approach to 
infinity and is equal to: 



(15) 
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6An(2n+ l)TO2i?^ 



> 0. 



For very large m, or small when the second term in the coefficient of the second 
derivatives in cqs. (12) and (14) can be neglected, numerical solution demonstrates that 
R would reach infinity in finite time in accordance with the results presented above (see 
Fig. 3, left panel). Nonzero g would terminate the unbounded rise of i?. To avoid too 
large deviation of R from the usual gravity coefficient g should be larger than or of the 
order of unity. In the right panel of Fig. 3 it is clearly seen, that for g = 1 the amplitude 
of oscillations remains constant whereas the average value of R increases with time. As 
follows from eq. (14), the frequency of small oscillations of y around yo = 1 + i^iTi in 
dimcnsionless time ri is 



(16) 



5 1+5^0^"+^ 
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Fig. 3. - Numerical solutions of eq. (14) for n = 3, ni — 0.01, y{Ti„) = 1 + niTin, y' {jin) = 0. 
Left panel: g = 0. Right panel: 5 = 1. 



It means that in physical time the frequency would be 

(17) --f(f) ; ^" 

In particular, for n = 5 and for a galactic gas cloud with Tq/Rq = 10^, the oscillation 
frequency would be 10^^ Hz « 10~^ eV. Higher density objects e.g. those with p = 
1 g/cm'^ would oscillate with much higher frequency, saturating bound (17), i.e. u ^ m. 
All kind of particles with masses smaller than m might be created by such oscillating 
field. 

3. Cosmological evolution and particle production in gravity 

In the present secton we study the cosmological evolution of the Universe in a theory 
with only an additional R^ term in the action, neglecting other terms which have been 
introduced to generate the accelerated expansion in the contemporary universe [13]. The 
impact of such terms is negligible in the limit of sufficiently large curvature, \R\ S> \Ro\, 
where Rq is the cosmological curvature at the present time. 

In other words, we study here the cosmological evolution of the early and not so early 
universe in the model with the action: 

(18) S ^ / d^xV-g (r-^)+ Sm . 
The modified Einstein equations for this theory read: 

(19) Rfj,i, — -gf_,i,R — ^— ^ [ Rfji, — -^RQiiu + 9ni^DaD°' — D^I?^ J R = -^^^^T^^^y . 



''PI 



Expressing the curvature scalar R through the Hubble parameter H = a/a as R 
—6H — 12H^, we get the time-time component of eq. (19): 



2H 2 3mj,iH' 
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where over-dots denote derivative with respect to physical time t. 
Taking trace of eq. (19) yields: 



(21) 



R + 3HR 



{ R 



In what follows, we study the cosmological evolution in the i?^-theory assuming rather 
general initial conditions for R and H and dominance of relativistic matter with the 
following equation for the matter content: 



(22) 



p + 4Hp = . 



It is convenient to rewrite the equations in terms of the dimensionless quantities 
T = Hot, h = H/Hq, r = R/Hq, y ~ Snp/ {3m^iHQ), and oj = m/iJo, where Hq is the 
value of the Hubble parameter at some initial time <o- Thus the following system of 
equations for dimensionless Hubble parameter is obtained: 



(23a) 
(23&) 



h" + 3hh' 



2h 



y 



2 h 

y' + 4ft.J/ = . 



First we assume that the deviations from General Relativity (GR) are small and 
expand h = l/(2r) + hi and y = l/(4r^) + j/i, assuming that hi/h <C 1 and yi/y <^ 1, 
and solve the linearized system of equations. 

The complete asymptotic solution for h has the form: 



(24) 



1 ci sin(wT + (/3) 
Mr) ^ + -^73 



The Hubble parameter oscillates around GR value, ^ 1/(2t) with rising amplitude, 
hi/ ho ~ T^^^, and for sufficiently large r the second term would start to dominate and 
the linear approximation would no longer hold. Using trancatcd Fourier expansion it is 
possible to obtain approximate analytical solutions of the full non-linear system in the 
high-frequency limit wr 3> 1. The same results are found numerically for the initial 
conditions /lo = 1 + Sh^ , h'^ — —2 + Sh'^ yo = 1 + Syo (see Fig. 4). 



0.50005 : |1 || I ,1 





Fig. 4. - Numerical solution of eqs. (23). Left panel: linear regime with 5ho = IQ'^, Sh'o = 0, 
yo = 1, uj = 10. Right panel: high frequency limit with 5ho — 1.5, Sh'o — 0, yo — 0, uj — 100. 
Initial conditions are different from GR, the central value hr = 0.6 is shifted from GR value 0.5. 
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Gravitational particle production may non-trivially affect the solutions of the above 
equations. Below we consider particle production by the external oscillating gravitational 
field and present the equation of motion for the evolution of R with the account of the 
back-reaction from particle production. This leads to an exponential damping of the 
oscillating part of R, while the non-oscillating "Friedmann" part remains practically 
undisturbed. The particle production influx into the cosmological plasma is estimated 
in the case of a massless, minimally-coupled to gravity scalar field with the action: 

(25) d^xV^g'^^d^clyd,^ . 

It terms of the conformally rcscalcd field, x = fl(*)<Aj ^-nd conformal time 77, such that 
adr] = dt, we can rewrite the equations of motion as: 



(26a) 


i?"4 


-2-i?' + r 
a 


(266) 


R = 


-6a" /a^ , 


(26c) 


X - 





1 



a 



c = o, 

We derive a closed equation for R taking the average value of the x-dependent quan- 
tum operators in the r.h.s. of eq. (26a) over vacuum, in presence of an external classical 
gravitational field R following the procedure described in ref. [14], where such equation 
was obtained in one-loop approximation. 

The dominant contribution of particle production is given by equation: 

(27) R + 3HR + m^R ^ /* dt' ^ . 
^ ' 12nm%Jt„ t-t' 

This equation is linear in R and naturally non-local in time since the impact of particle 
production depends upon all the history of the evolution of the system. 

Using again the procedure of truncated Fourier expansion including the back-reaction 
effects in the form of equation (27), we obtaine the decay rate: 

(28) Tr - 



This result is in agreement with ref. [15]. The characteristic decay time of the oscillating 
curvature is 

29 TR^ — = seconds. 

21 R rrf^ \ m J 

The contribution of the produced particles into the total cosmological energy density 
reaches its maximum value at approximately this time. 

The influx of energetic protons and antiprotons could have an impact on BBN. Thus 
this would cither allow to obtain some bounds on m or even to improve the agreement 
between the theoretical predictions for BBN and the measurements of the primordial 
abundances of light nuclei. 
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The oscillating curvature might be a source of dark matter in the form of heavy 
supersymmetric (SUSY) particles. Since the expected light SUSY particles have not yet 
been discovered at LHC, to some people supersymmetry somewhat lost its attractiveness. 
The contribution of the stable lightest SUSY particle into the cosmological energy is 
proportional to 

(30) VL ^ m%jjsY/mpi 

and for msuSY in the range 100 — 1000 GeV the cosmological fraction of these particles 
would be of order unity. It is exactly what is necessary for dark matter. However, it 
excludes thermally produced LSP's if they are much heavier. If LSP's came from the 
decay of R and their mass is larger than the "mass" of R, i.e. m, the LSP production 
could be sufficiently suppressed to make a reasonable contribution to dark matter. 

In contemporary astronomical objects oscillation frequency could vary from m down 
to very low frequency. The oscillations may produce radiation from high energy cosmic 
rays down to radio waves. 

^ ^ ^ 
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